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Abstract 



We study gauge mediation models of supersymmetry breaking with neutralino 
LSP. These models are naturally realized by embedding the usual four-dimensional 
gauge mediation models into a higher-dimensional spacetime such as M'^ x jZ^- 
We calculate the relic abundance of the neutralino LSP in these models and show 
. that there exist wide parameter regions where the neutralino LSP constitutes the 

dominant component of the cold dark matter. These regions evade constraints from 
collider experiments such as Higgs mass bounds and 6 — > 57, and also provide the 
value for the muon anomalous magnetic moment which is consistent with the SUSY 



5^ \ explanation of the deviation from the standard model prediction. 



1 Introduction 



Low-energy supersymmetry (SUSY) with dynamical SUSY breaking (DSB) is a very attrac- 
tive framework for explaining the large hierarchy between the electroweak and the Planck 
scales. One of the interesting features of this framework is that the lightest supersymmet- 
ric particle (LSP) is completely stable in models with R parity, which is well motivated 
by the stability of the proton. In conventional hidden-sector SUSY-breaking scenario, the 
LSP is believed to be the lightest neutralino, providing a good candidate for the cold dark 
matter in the universe [|l|. It is, indeed, suggested that the presence of a stable particle of 
weak-scale mass with electroweak interactions is a crucial ingredient for a natural solution 
to the cosmic coincidence problems 0. In fact, in view of its attractiveness, the possibility 
of neutralino dark matter has been extensively studied within the context of the hidden- 
sector scenario such as constrained version of the minimal SUSY standard model (MSSM) 

However, it is known that the hidden-sector scenario has no firm theoretical foundation. 
In this scenario, the Kahler potential is assumed to take a canonical form, but it can 
hardly be justified in supergravity. In supergravity, we expect the presence of the following 
non-renormalizable interactions in the Kahler potential: 

= ||ztz Qlg„ (1) 

which cannot be forbidden by any symmetry of the theories. Here, Z is the superfield 
responsible for the SUSY breaking and Qi represents generic standard-model quark and 
lepton superfields with z,j = 1, ■ ■ ■ ,3 being flavor indices; Mq — 2.4 x 10^^ GeV is the 
gravitational scale and rjij are constants of order one. These operators invalidate the as- 
sumption of the hidden-sector scenario, say the boundary condition of the constrained 
MSSM. Therefore, in the absence of any specific realization of the hidden-sector scenario, 
the mass spectrum used in most analyses of neutralino dark matter must be regarded as 
an ad hoc working hypothesis. [] Moreover, since there is no reason for the above operators 
to be flavor universal, they generically induce too much flavor-changing neutral currents 
(FCNC's) in the SUSY standard-model sector. In fact, suppressing FCNC is one of the 
most important issues in the SUSY standard model. 

Gauge-mediated SUSY breaking (GMSB) |rT| provides an elegant solution to this 

SUSY FCNC problem. In GMSB models, the squark and slepton masses are generated 

^ The special case that toq = ^1/2 = can be realized in five-dimensional theories with all the 
standard- model fields living in the bulk j^. 
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by the standard-model gauge interactions and are automatically flavor universal. On the 
other hand, in most GMSB models, the SUSY-breaking scale is supposed to be much lower 
than that in the conventional hidden-sector scenario. As a result, the gravitino mass m^j2 
becomes much smaller than the electroweak scale and the neutralino is no longer the 
LSP. The reason is again due to the non-renormalizable interactions Eq. (|I]). Since these 
operators would violate flavor universality, they generate flavor non-universal pieces of the 
soft SUSY-breaking masses for the squarks and sleptons of the order of the gravitino mass 
m^l2 — Fz/Mq- To preserve the success of GMSB models, these flavor non-universal pieces 
should be much smaller than the flavor universal pieces generated by the gauge mediation, 
which forces the gravitino to be much lighter than the electroweak scale. 

The question then is whether there are explicit models which solve FCNC problems 
and also accommodate neutralino dark matter. Notice that the above argument is entirely 
based on the presence of the non-renormalizable operators in Eq. (|l|), and it comes from 
the expectation in four-dimensional effective fleld theories that all operators consistent with 
the symmetries are present with coefficients of order one suppressed by some cut-off (grav- 
itational) scale. However, this naive expectation does not necessarily hold if fundamental 
theories are higher dimensional. Specifically, if Z and Qi fields are localized on different 
(3 + l)-dimensional branes in higher- dimensional spacetime, there is no direct contact in- 
teraction between Z and Qi fields and the fiavor non-universal operators in Eq. (|l|) are 



exponentially suppressed by the distance of two branes |]T2|, 0. Therefore, it is possible 



that the situations for the hidden-sector and GMSB models are drastically changed by em- 
bedding these models into higher dimensional theories at high energy scales. However, it 
turns out that the hidden sector still does not work, since necessary fiavor universal pieces 
for the squark and slepton masses are also suppressed in this case In the GMSB case, 

in contrast, the squark and slepton masses are generated by the standard- model gauge in- 
teractions. Thus, if we somehow manage the separation of Z and Qi, it merely says that 
the gravitino mass no longer has to be smaller than the electroweak scale. 

In fact, in a class of GMSB models [llO| , |15[| the sector responsible for DSB can be 



fully separated from the sector that feels standard-model gauge interactions; two sectors 
communicate with each other only through the U(l)m gauge interaction called messenger 
gauge interaction. Therefore, it is interesting to interpret these GMSB models as low-energy 



One possibility is to put the standard-model gauge fields in the bulk, and generate the squark and 
slepton masses through renormalization-group effects below the compactification scale . We here con- 
sider the case where the standard-model gauge and Higgs fields are also localized on the same brane as 
Oz's. 
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manifestations of the following brane- world scenario |16[. All fields in the DSB sector live 
on the DSB brane while the messenger and the standard-model fields are localized on our 
observable brane. The U(l)m gauge multiplet is put in the bulk, through which two sectors 
on different branes can communicate with each other. Then, the SUSY breaking on the DSB 
brane is transmitted to the observable brane only by the U(l)m gauge and gravitational 
interactions across the bulk.0 The non-renormalizable operators in Eq. (|lD are absent 
(exponentially suppressed) and flavor non-universal soft masses are not generated. 



Let us see how the above prescription works explicitly. It has been claimed |T8|, ^ that 
the brane separation produces the following no-scale type Kahler potential in low-energy 
four-dimensional effective field theories:^ 

K = -Slog (l - ^/o($obs, $1) - ^/d($dsb, $1)sb)) • (2) 

Here, $obs and $dsb denote superfields in the observable and DSB sectors, respectively. 
With the above Kahler potential, all soft SUSY-breaking masses and A terms in the ob- 



servable sector vanish at the tree level in the limit of zero cosmological constant [20, 13 



Then, the soft SUSY-breaking masses are generated only by the gauge mediation caused 



by the loop effects of the bulk U(l)m gauge interaction |jT6[. On the contrary, the term 
naturally arises at the tree level |^ from the Kahler potential if fo contains fo D HuHd, 



where Hu and Hd are chiral superfields of Higgs doublets. This mechanism |^ produces 
the SUSY-invariant mass fi of the order of the gravitino mass, i.e. fi c:^ fns/2- Therefore, if 
we choose the gravitino mass to be 771^/2 — 100 GeV — 1 TeV, we can correctly reproduce the 
electroweak symmetry breaking.^ Putting the U(l)m in the bulk, indeed, the transmission 
of the SUSY breaking becomes necessarily weak compared with purely four-dimensional 
GMSB models. Thus, we can naturally obtain the gravitino mass — 100 GeV — 1 TeV 
with moderate size of extra dimensions. In the five- dimensional model of Ref. for 
instance, the compactification length L is given by L ~ (2 x 10^^ GeV)~^, which is close 
to the value obtained in the scenario of Ref. []T^|. This provides a simple solution to the 



fi problem in GMSB models with DSB and observable sectors geometrically separated in 



higher-dimensional spacetime 17, 16 



The consequence of the above brane-world GMSB models is that the gravitino has a 



Another mechanism of transmitting SUSY breaking between two branes is discussed in Ref. [ p7[ . 

^ The no-scale supergravity adopts a specific form fo = Z + . We here assume the Kahler 
potential for the Z field to be of the form fD = Z^Z + ---, where the ellipsis denotes higher order terms. 

^ With these gravitino masses the anomaly mediation generates too small SUSY-breaking masses 

in the observable sector. 
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mass of order 100 GeV — 1 TeV, although the mass spectrum of the gauginos, squarks 
and sleptons is the same with that of the usual GMSB models [0, |1^. Thus, the lightest 
neutralino is most likely the LSP and can behave as a cold dark matter in the universe. 
In this paper, we calculate the relic abundance of the neutralino LSP in these brane- 
world GMSB models and show that the neutralino can actually be the dark matter in a 
wide range of the parameter space. Since GMSB models are highly predictive, the relic 
abundance is calculated in terms of a few parameters: one parameter smaller than in 
the case of constrained MSSM. Requiring that is in the cosmologically favored region, 
0.1 ^ ^ 0.3 {h is the present day Hubble parameter in units of 100 kms~^Mpc~^), 
we obtain the constraint on the parameter space of the models and hence masses for the 
superparticles. We also calculate the lightest Higgs boson mass and the constraint from 
b s'-f process to identify the phenomenologically allowed region of the parameter space. 
The resulting region is remarkably consistent with the SUSY explanation of the recently 
reported 2.6a deviation p3| of the muon anomalous magnetic moment from the standard- 
model value. We assume the conservation of R parity throughout the paper. 



2 Gauge-Mediated SUSY Breaking Model 

The mass spectrum of the superparticles has a great influence on the estimation of the 
LSP relic abundance. Since brane-world GMSB models have the same mass spectrum as 
ordinary four-dimensional GMSB models except for the gravitino, we begin with reviewing 
the soft masses for the gauginos and sfermions in GMSB models. 

In GMSB models, the superparticle mass spectrum does not depend on the detail of the 
DSB sector. Thus, it is sufficient to simply consider the messenger sector which consists 
of N pairs of vector-like messenger superfields and (ji {i = 1, ■ ■ ■ , N). To preserve the 
gauge coupling unification, is supposed to form a complete grand unified theory (GUT) 
multiplet. We assume, in this paper, that and qi transform as 5 and 5* representation 
under the SU(5)gut- Then, the superpotential for the messenger sector is written as 

N _ _ 

w = J2 {>^tsdA + Kskh) , (3) 

i=l 

where {di, U) G qi and (rfj, U) G qi. The messenger quark multiplets d and d transform as 
the right-handed down quark and its antiparticle under the standard-model gauge group, 
respectively, and the messenger lepton multiplets / and / as the left-handed doublet lepton 
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and its antiparticle, respectively. The singlet field is a spurion superfield parameterizing 
the SUSY-breaking effect, and assumed to have vacuum expectation values (VEV's) in its 
lowest and highest components as (S) = M + 

The SUSY breaking effect in the messenger sector is transmitted to the superpartners of 
the standard-model particles through the standard-model gauge interactions. The gaugino 
masses and the sfermion squared masses are generated at the messenger scale, M, by one- 
and two-loop diagrams, respectively. In general, various coupling constants, \f and in 
the superpotential Eq. are not equal at the messenger scale; for instance, the couplings 
of the messenger quarks and leptons are not the same at the scale M, X'^ ^ X\ even if we 
assume they are equal at the GUT scale. However, the effect of having different couplings 
on the superparticle mass spectrum does not appear at the leading order in F/M"^, so that 
we here take all the coupling constants in Eq. (|^) to be equal for simplicity and absorb it 
into the definition of M and F. Then, the gaugino masses Ma and the sfermion squared 
masses are given by p^, p5| : 

F 



m 



f 



AttM 
M 



M2 



2N 



M2 



(4) 
(5) 



where a = 1, 2, 3 represents the standard-model gauge groups and is the quadratic 
Casimir coefficient for the representation each sfermion belongs to. Here, the functions Q 
and J-' are defined by 



g{x) 



— \{l + x) ln(l +x) + {l-x) ln(l - x)] 
x^ 

[ln(l + x) -2Li 



x" 



2x 



+ (x 



-X] 



(6) 
(7) 



These functions have properties that Q{x) ^ ~ 1 when x ^ 1. 

To obtain the mass spectrum of superparticles, we have to evolve soft masses given in 
Eqs. (^, ^ using renormalization-group equations (RGE). Together with the contribution 
from the electroweak symmetry breaking, the masses for all the superparticles are deter- 
mined. The Higgs sector contains two more parameters yU and ^B: the SUSY-invariant 
{W = fi HuHd) and SUSY-breaking {C = fxB huhd) mass terms for the Higgs doublets. In 
brane-world GMSB models, both /j, and B are generated of the order of the weak scale, in 
contrast to the minimal model discussed in Ref. [EgI where B = is assumed at the messen- 



ger scale. The electroweak symmetry breaking condition relates these two parameters to v 
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and tan/3 up to the sign of /i, where v = \J {huY + {hdY — 175 GeV and tan/? = {hu) / {ha). 
Therefore, we end up with the following set of free parameters in our analyses: N, F/M, 
M, tan/5 and sgn(/i). 

With the above GMSB mass spectrum, the correct electroweak symmetry breaking 
requires rather large value of the /i parameter, /i^ 2M2. This is because in GMSB models 
the colored particles are relatively heavy and, as a result, the Higgs-boson mass squared 
receives large negative contribution from the top squark through the top Yukawa coupling. 
This fact has some important consequences in the present brane-world GMSB scenario. 
First, since /i is of the order of the gravitino mass, the gravitino tends to be heavier than 
the lightest superpartner of the standard model. Thus, the gravitino is not the LSP in 
contrast to the usual four- dimensional GMSB models. According to Eqs. then, 
the lightest neutralino x o^' the right-handed stau f/j can be the LSP. In this paper, we 
concentrate on the case where x is the LSP and is the next to LSP (NLSP), since f^i 
LSP leads to the serious problem of charged dark matter. Indeed, the LSP is x most of 
the parameter space, especially when = 1. 

Another important consequence of /i ^ M2 is that the lightest neutralino is almost 
purely composed of the bino. This leads to a significant simplification in understanding the 
neutralino relic abundance; for instance, annihilation into Zh is strongly suppressed due to 
the smallness of the Higgsino component in x- In the next section, we list the processes 
relevant for determining the relic abundance of the lightest neutralino aiid calculate 
various quantities in the present models. 

3 Relic LSP Abundance 

The LSP is stable in i?-parity preserving models, and hence its number density can de- 
crease only through annihilation processes. In an expanding universe, the pair-annihilation 
"freezes out" when the expansion rate of the universe exceeds the interaction rate of the 
annihilation. After the freeze out, the number density of the LSP per comoving volume is 
constant, so that some amount of relic LSP is left in the present universe. Neutralino LSP 
abundance has been well studied in the context of conventional hidden-sector SUSY break- 
ing models p, ||, ^ H, 0. It is easily estimated once we could determine the relevant cross 
sections for the annihilation of the LSP. We first briefly review the procedure of calculating 
the relic abundance. 

The time evolution of the LSP number density is described by the Boltzmann equation. 
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Since in GMSB models the bino and the right-handed stau are degenerate in some parameter 
region, we have to take into account the coannihilation effect [B^, EHl. The Boltzmann 



equation with the coannihilation effect can be written as an equation for n = J^i^i^ where 
rii are the number densities of the species i, and i represents the LSP neutralino x and the 
right-handed charged sleptons tr, f^, fiR, fi*j^, cr and e\. The equation takes the following 
form 0: 

dn 
It 

where H is the Hubble parameter and 



— = -3Hn - {(Tcsv) 



n 



[n 



eq\2 



^eq ^eq " 



(9) 



Here, n'^'^ (^i'^) is the equilibrium value of n (nj) and v is the relative velocity of the particles 
i and j. The bracket denotes the thermal average and aij is the total annihilation cross 
section of i + j —>■ X + X': 



X,X' 



X + X') 



(10) 



where X and X' represent possible standard-model particles. 

There are 7x7 cross sections {<Ji/s) in Eq. (H), but most of them are not indepen- 
dent. The independent cross sections are listed in Table [l|, where we have shown only 
relevant final states which are kinematically accessible and have non-negligible cross sec- 
tions. Furthermore, there are some simplifications coming from the GMSB mass spectrum. 
For XX ~^ f f 1 foi' example, the cross section is dominated by IrIr final states since the 
right-handed sleptons are much lighter than the other sfermions and have the largest value 
of the hypercharge. 

To estimate the relic LSP abundance, we need {a^sv) at the freeze-out temperature T/ 
of the LSP. Since typically Tf ~ where is the LSP mass, the expansion of {a^^v) 

in terms of T /my. (partial-wave expansion) is relevant. Thermally-averaged cross section 
{aijv) for the process i + j ^ k + I is given by : 



m^m^ 



3{mi + mj 
2mimj 



-T 



*{mi+mjy+3{mi+mj)T 




where 



wis) = 



{2'Ky5\pi+p,-pu-pi)\r\ 



(12) 
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Initial state 


Final states 


XX 


// 


I'rX 


/S, I'z, rh 


li li* 




nJ' Ft 


rP 


ly^Ri^^j) 





Table 1: Annihilation cross sections; z,j = r, /i, e. 



Here, s = {pi + PjY is the Mandelstam variable, and |Tp is the transition matrix element 
squared summed over final state spins and averaged over initial state spins. Terms of order 
(T/m^)^ and (T/m^Y are called s-wave and p-wave components, respectively. 

The s-wave component of the thermally-averaged neutralino annihilation cross section 
(cr^^v) is suppressed by tiny final-state fermion masses, so that the p-wave part is the dom- 
inant piece. While neutralino annihilation cross section is p-wave suppressed, those for the 
sleptons, {aj.f.v), have s-wave components as dominant pieces. Therefore, if the equilib- 
rium number densities for the sleptons are not much smaller than that for the neutralino, 
the slepton annihilation processes can significantly reduce the relic LSP abundance (see 
Eq. ®). This coannihilation process is effective when the slepton masses are degenerate 
with the neutralino mass within ~ 20%. In the GMSB spectrum, it happens when tan j3 is 
large and/or is greater than 1. 

The partial-wave expansion of the thermally-averaged cross section does not give a 
good approximation when the initial momentum is near the s-channel pole or the final- 
state threshold |2^. They could occur at Z, h, A, H poles and WW, ZZ, Zh, ti thresholds 
in the neutralino annihilation. However, with the present GMSB mass spectrum, these cases 
simply do not happen or their effects are strongly suppressed. First, the Z, h poles can be 
hit when m<^ ~ mz/2,mh/2, but these regions are already being excluded by the chargino 
search at LEP2. The situations rriy. ~ mA/2,mH/'2 also do not occur in the parameter 
region we are considering. As for the threshold effect, WW and ZZ ones are small due 
to /i ^ Ml] Zh one is also negligible due to the smallness of the Higgsino component in 
X, and tt one is strongly suppressed by the large masses for the top squarks exchanged in 
t-channel. 

With these understandings, we can calculate the relic abundance Vt^h"^. In the actual cal- 
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culation, we have used the computer program neutdriver coded by Jungman, Kamionkowski 
and Griest |p, which contains all the (co) annihilation cross sections calculated by Drees 
and Nojiri In Fig. |l|, we have shown a cosmologically favored region, 0.1 < f^^^^ ^ ^-^ 
(light shaded regions), on M-F/M plane in the case of (A^, tan/3) = (1, 10). We scanned 
the region 10^ GeV < F/M < 2 X 10^ GeV, which corresponds to the soft SUSY-breaking 
masses ~ 100 GeV - 1 TeV. The sign of fi is taken to be positive in the standard notation 
(the one in which the constraint from 6 ^ 57 process is weaker). The region extends from 
upper left to lower right directions. This is because the relic abundance is almost com- 
pletely determined by the mass of the right-handed stau, which is monotonically increasing 
with M with a fixed value of F/M due to renormalization group effects. We have also 
drawn the contours of the lightest Higgs boson mass (solid lines) and the lightest chargino 
mass (dashed lines) in GeV, which are calculated using neutdriver. We find that some of 
the parameter region satisfies constraints from the lower bounds on the Higgs boson mass 
rriho ^ 113.5 GeV and the lightest chargino mass m^± ^ 150 GeV. In particular, we find 
that smaller messenger scale, M ^ 10^ GeV, is favored. 

In Fig. D, we have shown the constraint from b 57, 2.3 x 10~^ < Br{b 57) < 
4.1 X 10~^ ||2^, in the same M-F/M plane as Fig. The dark shaded region indicates the 
excluded region. We see that most of the parameter region which satisfies the constraints 
from the Higgs and chargino masses also satisfies that from b — > s'j. We have also drawn the 
contour of the SUSY contribution to the muon anomalous magnetic moment, a^, in units of 



10 . The SUSY contributions to were discussed, for example, in Refs. |30|, ^ It 



is interesting that the value is consistent with that required to explain the recently claimed 
2.6a deviation of between the observed and the standard-model values. 

In Figs. 1^ and we have plotted various quantities in the case of (A^, tan/3) = (1, 50). 
In this case, the coannihilation effect is important so that the cosmologically favored region 
(light shaded region) is shifted to larger values of F/M in lower M region (M ~ F/M). We 
find that the constraint from 6 — 57 is more stringent than the tan/3 = 10 case, but the 
region 10"^ GeV still satisfies the constraints. This region gives the value of a^|susY 
within la and 2a level for the SUSY explanation of the observed deviation. 

We have also done the same analyses in the case of = 2, which are plotted in 
Figs. 1^, and |[ The black region represents the region where the stau is the lightest 
SUSY particle. In the case of tan/5 = 10 there is a consistent region which reproduces 
cosmologically interesting abundance of the LSP, but we do not find such a region in the 
case of tan /3 = 50. 
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4 Conclusions and Discussion 



In this paper we have investigated gauge mediation models with neutrahno LSP. This type 
of models is naturally realized in the extra-dimensional setup |T^. Due to the geometrical 
separation between the SUSY-breaking and observable sectors, flavor non-universal con- 
tributions to the squark and slepton masses are exponentially suppressed, which makes it 
possible to solve the fi problem by the mechanism of Ref. [^. The consequence is that 
the gravitino mass becomes of the order of the weak scale, while all the other superparticle 
masses are the same as those in usual GMSB models. Therefore, the lightest neutralino 
(mostly bino) could provide the cold dark matter of the universe. 

We have calculated the relic abundance of the LSP neutralino in this brane- world GMSB 
model. We found that there is a wide parameter region where cosmologically favored relic 
densities, 0.1 < Q^.h'^ < 0.3, are obtained and, at the same time, the lower bounds on 
the Higgs and the lightest chargino masses from collider experiments are evaded. We 
have also found that the experimental constraint on the 6 — 57 process is satisfied in 
this cosmologically favored parameter region, especially when the messenger scale M is 
low. Interestingly, this parameter region also gives the SUSY contribution to the muon 
anomalous magnetic moment, a^lsusYj with the values required to explain the recently 
reported 2.6a deviation of p3[ from the standard-model value. These are due to the 
fact that in the GMSB spectrum the colored superparticles are relatively heavy while non- 
colored ones are light; the constraint from 6 — > 57 is easily evaded since squarks are heavy, 
while we obtain relatively large values of a^|susY since sleptons are light. Although we 
have limited ourselves to the case of the minimal messenger sector in this paper, the above 
features are generic to GMSB models, for example, ones with more general messenger 
sectors [^, ^ |3^. We leave the study of the LSP relic abundance under these more 
general superparticle spectrums for future work. 

It is interesting to compare the present situation with other SUSY-breaking models that 



solve the SUSY flavor problem. In anomaly mediated SUSY-breaking scenario [|1^, |22[ , the 
LSP is generically wino-like neutralino, and its annihilation cross section is too large to 
obtain cosmologically favored thermal relic density of the LSP, so that we need to consider 



some other source of the LSP, such as the decay of the moduli field In the decoupling 
scenario MTl and the focus-point scenario it would be difficult to accommodate suffi- 



ciently large values of the muon anomalous magnetic moment since the smuons are assumed 



to be heavy in these scenarios. The gaugino mediation scenario [14] is a promising candi- 
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date theory, but the parameter region of the neutrahno LSP is not necessarily very large, 
especially when the Higgs fields are localized on the observable brane. Thus, we conclude 
that the GMSB models with neutralino dark matter provide one of the promising SUSY- 
breaking scenarios which are in accordance with current phenomenological situations. It 
would be interesting to study further the relic abundance of the LSP neutralino and its 
detections in this context. 
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Figure 1: GMSB parameter space (M, F/M) corresponding to 0.1 < Vt^h^ < 0.3 (the light 
shaded region) for tan/3 = 10, sgn/x = +1. The dashed hnes 100, 200, 300, 400, 500 are the 
contours for the hghtest chargino mass, and the sohd hnes 100, 110, 120 are the contours 
for the Higgs boson mass (in unit of GeV). The dark shaded region in the upper left corner 
is the stau LSP region. 
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Figure 2: GMSB parameter space (M, F/M) corresponding to 0.1 < Vl^h"^ < 0.3 (the light 
shaded region) for tan/5 = 10, sgn/i = +1. The dotted hnes 11, 27, 43, 59, 75 are the 
contours of the — 2(T, — Icr, central, +1(T, and +2(T values of the SUSY contribution to the 
muon ((7 — 2), a^, in unit of 10~^°. The dark shaded region in the bottom is excluded by 
h 57. 
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Figure 3: The same figure as Fig. |T] for tan/? = 50. 
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Figure 4: The same figure as Fig. |^ for tan [5 = 50. Tfie narrow strip in tfie dark sliaded 
region (excluded hy b —>■ 57) is an accidentally allowed region. 
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Figure 5: The same figure as Fig. |T] for N = 2. The dark sfiaded region in tlie upper left 
(stau LSP region) is larger than that for = 1. 
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Figure 6: The same figure as Fig. ^ for N = 2. The dark shaded region in the upper left 
(stau LSP region) is larger than that for = 1. 
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Figure 7: The same figure as Fig. ^ for N = 2. In tliis case tlie region for 0.1 ^ i^^/i^ ^ 0.3 
is absent in the range of the parameters exhibited. 
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Figure 8: The same figure as Fig. ^ for N = 2. In tfiis case tfie region for 0.1 ^ VL-^h^ ^ 0.3 
is absent in tlie range of the parameters exhibited. 
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